Theory of optical conductivity for dilute Gai_xMnxAs 
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We construct a semi-microscopic theory, to describe the optical conductivity of Gai_xMnxAs in 
the dilute limit, a; ~ 1 %. We construct an effective Hamiltonian that captures inside-impurity 
band optical transitions as well as transitions between the valence band and the impurity band. All 
parameters of the Ifamiltonian are computed from microscopic variational calculations. We find a 
metal-insulator transition within the impurity band in the concentration range, x ~ 0.2 — 0.3 % for 
uncompensated and a; ~ 1 — 3 % for compensated samples, in good agreement with the experiments. 
We find an optical mass mopt ~ rrie, which is almost independent of the impurity concentration 
excepting in the vicinity of the metal-insulator transition, where it reaches values as large as mopt ~ 
10 me. We also reproduce a mid-infrared peak at hui ~ 200 meV, which redshifts upon doping, in 
quantitative agreement with the experiments. 

PACS numbers: 75.50.Pp, 78.66.Fd, 78.66.Fd, 78.20.-e, 78.67.-n. 
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I. INTRODUCTION 

The diluted magnetic semiconductor (DMS) 
Gai_a;Mn2;As has emerged as one of the most promising 
materials due to its potential applications in spin-based 
technology. This material is, however, equally interest- 
ing from the point of view of fundamental research for 
its unique properties and the rich physics it displays; 
Gai-xMuxAs is a strongly disordered ferromagnet, 
where the interplay of ferromagnetism, localization, 
magnetic fluctuations and the presence of strong 
spin-orbit coupling lead to many interesting proper- 
ties such as a strong magnetoresistancefijS resistivity 
anomalies j^i^i^ or the presence of a large anomalous 
Hall effect^ 

Although Gai_xMnxAs has been the subject of intense 
theoretical and experimental investigation, even its most 
basic properties are still debated. One of these funda- 
mental and unresolved issues is the existence or non- 
existence of an impurity band in this material. There is a 
general consensus that for very small magnetic concentra- 
tions, X < 1%, Gai_xMnxAs is described in terms of an 
impurity band. For higher concentrations, however, there 
is no general consensus yet: On one hand, essentially all 
spectroscopic measurements such as ARPES^' STM,^ op- 
tical conductivity data^ i^°'^^i^^ and ellipsometryH seem 
to favor the presence of an impurity band even up to 
moderate concentrations (a; « 3 — 5%), and even high 
quality samples with a high Curie temperature and a 
clearly metallic behavior have surprisingly small values 
of kpl ^ Ij^ii^ On the other hand, many properties of 
these materials can be even quantitatively understood in 
terms of a disordered valence band picture, and, in fact, 
from a theoretical point of view, it is hard to understand 
how an impurity band could survive up to concentrations 
as high as a; « 5%fi^ 

In the present paper we shall not attempt to resolve 
this discussion, rather, we would like to focus exclusively 



on the optical conductivity. Typical optical conductiv- 
ity measurements in Gai_xMnxAs display two rather re- 
markable features: (i) a mid infrared resonance peak at 
approximately 200 meV that redshifts with increasing 
hole concentration, p; (ii) a large optical mass of the or- 
der TTiopt ~ 0.7 1.5me, implying a mobility which is 
orders of magnitude smaller than the one of GaAs doped 
with similar concentrations of non-magnetic impuritiesJ^ 

In this paper we aim at developing a semi-microscopic 
theory for the optical conductivity of Gai_xMnxAs in the 
very dilute limit, x ^ 1%. There are several theoretical 
studies of the optical conductivity of Gai_xMnxAs by 
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Maybe the most realistic calculations 
have been carried out in Refs.[l^ and however, 
the results of these works are somewhat conflicting: In 
Ref. [3| a disordered valence band described in terms of 
the Luttinger model has been studied. In these calcula- 
tions, rather surprisingly, the mid-infrared peak appeared 
even in a simplified parabolic model, where a completely 
isotropical valence band mass was assumed and the spin- 
orbit split valence band is completely ignored. Therefore 
the authors concluded that the mid-infrared peak is not 
due to transitions to the spin-orbit split valence band, 
rather, they interpreted it as a result of strong localiza- 
tion and interference effects. In another, very thorough 
calculation^^, realistic tight binding models have been 
used, but the results were conflicting in that features as- 
sociated with an impurity band have or have not ap- 
peared depending on the particular tight binding scheme 
applied. 

Unfortunately, none of these previously applied meth- 
ods is appropriate to study accurately the dilute limit, 
X ^ 1%. This small concentration limit is of particular 
interest, since the metal- insulator transition takes place 
at concentrations as low as x « 0.5 and further- 

more, it is in this limit where features associated with 
the impurity band should be present. In the present pa- 
per, we shall follow the lines of Ref. [l^ and develop 
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a semi-microscopical approach to capture the physics of 
this dilute limit. At the microscopic level, we describe the 
valence band holes as spin F ^ 3/2 carriers, which in- 
teract with the external electromagnetic vector potential 
and the Mn ions, as described by the following Hamilto- 
nian: 

i i,m 

+ ^ J ,5(r, - R„0 F, • S„ . (1) 

Here and Pi {i = 1, . . . , Nh) denote the valence hole 
coordinates and momenta, and the Mn atoms are located 
at random positions, R„i (i — 1, . . . , -/Vmh)- In the first 
term we assume an isotropic effective mass toq — 0.56 rrie. 
This is a somewhat uncontrolled approximation. How- 
ever, according to the results of Ref. [lB|, inclusion of 
a more realistic band structure influences only slightly 
the optical spectrum for frequencies fioj < SOOmeV, we 
shall therefore ignore it here. The second term describes 
the random potential created by the Mn atoms, and it 
incorporates the Coulomb attraction as well as the cen- 
tral cell correction.™ Finally, the last term describes the 
local exchange interaction between the manganese mag- 
netic moments, Sm and the spin of hole i, F^. The pa- 
rameters of the 'microscopic' Hamiltonian Eq. ([1]) are all 
well-known. 

In our work we use the Hamiltonian ([1]) as a start- 
ing point, and derive a microscopic effective Hamilto- 
nian from it that captures the impurity band physics 
in the very dilute limit, but also accounts for impu- 
rity band-valence band transitions (for details see Sec- 
tion |Tl|. This effective Hamiltonian can then be used to 
compute the optical conductivity using field theoretical 
methods, as detailed in Section IIIII With this 'multi- 
scale approach', we are then able to determine the op- 
tical properties of Gai^^MuxAs accurately in the small 
concentration regime, while keeping « 200. 

One of the key elements of our approach is the as- 
sumption that, for these small concentrations, the Fermi 
level lies inside the impurity band. This is indeed evi- 
denced by transport measurements;^ which clearly show 
that the activation energy (associated with localized im- 
purity states to valence band transitions) remains fi- 
nite as one approaches the metal insulator transition. 
This implies that the metal-insulator transition takes 
place within the impurity band, and one must be able 
to capture the properties of small concentration metal- 
lic Gai-xMuxAs samples using our impurity band-based 
method. 

Indeed, our results also support this picture; a detailed 
analysis of the hole states reveals a metal-insulator tran- 
sition at about a; « 0.2 0.3 %, where delocalized states 
appear in the middle of the impurity band. As shown 
in Fig. [1] this is approximately the critical concentration 
for uncompensated samples. The Anderson transition 
thus takes place inside the impurity band, in agreement 
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FIG. 1: Density of states on for x = 0.2% and x = 0.3%. 
Both the impurity band and valence band contributions are 
shown. At these concentrations the impurity band is clearly 
separated from the valence band. The gray area indicates the 
regime of localized states, while delocalized states are shown 
as white regions. For these calculations we used a hole fraction 
/ — 0.5. The system size is measured in the size of the unit 
cell, ao. 



with the conclusions of Ref. and also in good agree- 
ment with activation energy values observed for sam- 
ples grown at high temperature . ^^'^^ Typically, however, 
Gai-xMuxAs samples grown under non-equilibrium con- 
ditions are compensated even after annealing: in addition 
to the substitutional Mn ions (of concentration xs) there 
is also a finite concentration xi of interstitial Mn ions, 
which behave as double donors, and are also believed to 
make a fraction of the Mn ions inactive by simply bind- 
ing to them.^^ As a result, the effective concentration of 
"active" Mn ions is reduced to x = Xctt ~ xs ^ xj and 
the concentration of holes is also suppressed compared 
to XcS by the hole fraction / « (xs — 2xi)/{xs — xi). 
Remarkably, even if only 20% of the Mn ions goes to 
interstitial positions, that reduces the effective concen- 
tration to 60% percent of the total Mn concentration 
and amounts in a hole fraction / w 0.66%. Unannealed 
samples tend to have even larger interstitial Mn concen- 
trations and thus much smaller hole fractions f^^i^liS^ As 
a result, depending on the precise annealing protocol, 
ferromagnetic Gai_xMnxAs samples (grown under non- 
equilibrium conditions) tend to show the phase transi- 
tion at higher concentrations. According to our calcu- 
lations, for xz/xtotai = 30% (/ = 0.3) the metal insu- 
lator transition takes place at about xs ~ 3.2 %, while 
for 2;//a;totai = 23% (/ = 0.55) the transition occurs at 
Xs ~ 0.7 %. These values are consistent with the experi- 
mental datai^'^Si^S'^i^ In the rest of the paper, we shall 
not make a difference between substitutional and inter- 
stitial Mn ions, rather, by the symbol x we shall refer 
to the concentration of active / effective Mn ions, which 
participate in the formation of the ferromagnetic state. 
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FIG. 2: Comparison between our theory and experimental 
data. The optical conductivity data presented in the inset 
were extracted from Ref. II 111. 



and / will be the corresponding hole fraction. 

We also find that the optical mass increases to very 
large values Wopt ~ 10 rrie at the critical concentration, 
where the DC conductance vanishes. However, apart 
from this, the optical spectrum changes rather contin- 
uously. For larger Mn concentrations, a; > 2 — 3 %, we 
obtain an optical mass mopt ~ Trie, which is almost inde- 
pendent of the hole concentration p, and is in surprisingly 
good agreement with the experimental valueSfii 

Our approach is designed to work in the limit of small 
concentrations, and it should break down for large con- 
centrations. Where exactly this breakdown takes place, is 
not quite clear. Our calculations show that the impurity 
band is not merged with the valence band for concen- 
trations as large as a; ~ 3%. Furthermore, the impurity 
band picture may make sense for even larger concentra- 
tions if the Fermi level is inside the gap (as indicated 
by ARPES measurements),^ where states are expected 
to exhibit a stronger impurity state character. We are 
not able to convincingly determine the upper limit of 
our approach. However, if we blindly extend it to the 
regime x > 3%, without justifying their applicability, 
rather surprisingly, our calculations qualitatively as well 
as quantitatively reproduce all important features of the 
experiments. 

One of the typical optical conductivity results is dis- 
played in Fig. [21 where, for comparison, we also show 
some typical experimental data.^^ In our calculations the 
Drude contribution originates entirely from carriers re- 
siding within the impurity band, and the 'Drude peak' 
appears as a plateau at smaller frequencies, just as in 
the optical conductivity experiments. The mid-infrared 
peak, on the other hand, is due to transitions from the 
impurity to the valence band. Notice that the position 
as well as the overall size of the signal are quantitatively 
reproduced. 

Having the Fermi level within the impurity band natu- 
rally explains the redshift of this resonance with increas- 
ing hole concentration. In Fig. Owe show the computed 
peak position as a function of the effective optical spec- 
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FIG. 3: Mid-infrared peak position extracted from the opti- 
cal conductivity curves vs. the effective optical spectral weight, 
A^eff (defined through Eq. (f28)l '). Open symbols represent the- 
oretical results while the solid dots are the experimental data 
extracted from Ref. [llj. Calculations for different hole frac- 
tions / and concentrations x {x — 1 5%) fall onto a single 
curve. 



tral weight, A'off (defined through Eq. (UHl)) for a variety 
of concentrations x and hole fractions, / = Nh/NMn- 
Our calculations show excellent agreement with the ex- 
perimental results of Ref. Ull . 

The paper is organized as follows. First, in Section Ull 
we outline the calculations that lead from the microscopic 
Hamiltonian Eq. ^ to the effective Hamiltonians used. 
Some of the details of these rather technical variational 
calculations are given in appendices El and |B1 Sec. IIIII 
provides the basic expressions for the optical conductiv- 
ity within a mean field approach, while the results of our 
numerical calculations are given in Sec. IIVI In Sec. IVlwe 
show how to incorporate the effects of magnetic fluctua- 
tions, and finally, we conclude in Sec. IVII 



II. THEORETICAL FRAMEWORK 



With current day computer technology, it is essentially 
impossible to treat the Hamiltonian of Eq. ^ accu- 
rately for small Mn concentrations, x ^ 1%. We there- 
fore describe this regime in terms of an effective Hamil- 
tonian, which we derive from Eq. (jl]) , and which consists 
of an impurity band part and a valence band part. The 
first part of the present Section focusses on this map- 
ping. However, to compute the optical conductivity, we 
also need to determine how the electromagnetic field in 
Eq. ^ couples to impurity states within the effective 
Hamiltonian. This is discussed in the second part of the 
Section. In the main body of the text we only outline 
and summarize the most important steps. Details on 
these simple but rather lengthy calculations are given in 
Appendices [XllBl and are provided as supplementary in- 
formation. 
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A. Effective Hamiltonian 

1. Impurity hand Hamiltonian 

Since for very small concentrations the Fermi energy 
is in the impurity band, we first need to describe the 
impurity states. In the small concentration limit, hole 
wave functions are localized at the Mn sites, providing 
a strong and attractive potential for them. Therefore 
the impurity band can be described using a tight binding 
Hamiltonian of the form>2S,)2^ 



Here creates a bound hole at Mn position 
spin state \F = 3/2, = a), and tij denotes the hopping 
between Mn sites Ri and Rj. Since the hole mass is 
isotropic in our microscopic Hamiltonian, the hopping is 
independent of the hole spin a.^'^ The on-site energy Ei 
contains the binding energy of the hole, Eq = —110 meV, 
but it also accounts for the Coulomb and kinetic energy 
shifts generated by neighboring Mn ions. Finally, the last 
term describes the local antiferromagnetic exchange with 
the core Mn spin Si. 

The parameters appearing in Eq. ([2]) can be deter- 
mined from experiments and from the microscopic Hamil- 
tonian, Eq. Il]) . The coupling G is known from EPR ex- 
periments to be G ~ 5meV for a single Mn impurityj ^'^i'^^ 
However, the hopping parameters and the on site energy, 
Ei need to be computed from Eq. ^ . Similar to Ref . [23 . 
we determined them from a variational calculation of the 
molecular orbitals for an Mn2 'molecule', described by 
the Hamiltonian: 

1 1 



(2) 



m a 



Ho 



7, 



eri 



er2 



+ VcM) + Vcc{r2) . (3) 



where we used atomic units, ri 2 = |r — Ri 2I, e = 12.65 is 
the dielectric constant of GaAs, and 7 — irie/mo — 1.782. 
The so-called central cell correction Vcc accounts for the 
local interaction at the Mn core, and is given hy^ 

V{r) = Vcc{r) = -Vo e'^'''^ . 



The parameters Vq and ro must be chosen to repro- 
duce the experimentally observed impurity state at E^ k, 
llOmeV. In our calculations we have used Vb = \.'oeV 
and ro — 2A, but our results do not depend on this par- 
ticular choice as long as rp is sufficiently small. Details 
of this variational calculation are presented in Appendix 
iB]) . In these calculations we neglected the coupling G, 
since it is much smaller than the binding energy of the 
holes. The final result is that (for G = 0) the low ly- 
ing spectrum of two Mn ions at a separation R can be 
described by the following simple effective Hamiltonian, 



H = -t{R)J2icLc2c. + h.c.) + E{R) 



4=1,2 



(5) 
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FIG. 4: Top: Variationally computed hopping parameter t 
as function of the relative distance between the Mn sites, R. 
Bottom: The energy shift AE{R). The inset shows the cor- 
rection AE{R) to the classical Coulomb shift. 



The parameters t{R) and AE{R) = E{R) — Eo are shown 
in Fig. m As the inset of the lower panel of Fig. IHshows, 
much of the energy shift AE{R) originates from a simple 
long-ranged Coulomb shift due to the Coulomb poten- 
tial of the neighboring Mn site, AEcouiomb{R) = — 1/ei?. 
The remaining kinetic shift, AE{R) = AE{R) + 1/ei? is 
relatively large for small separations, but vanishes expo- 
nentially for large values of R. 

Having determined the spectrum and the effective 
Hamiltonian of the two-Mn ion complex, we can now 
express the parameters of the effective Hamiltonian, 
Eq. ([2]). The hopping parameter between sites i and j 
(4) depends just on their separation, Rij = [Ri — Rj|, and 



is simply tij = t{Rij), the hopping obtained from the 
two-Mn problem, while the on-site energies in Eq. ^ 
are given as 



E, ^ Eo + AE, , 

AE, = Y^AE{R,, 

3^1 



(6) 



The sum in this equation should be evaluated carefully: 
As we discussed above, for large separations, the shifts 
AEiRi — Rj) are dominated by the long-range Coulomb 
contributions, and therefore the sum in Eq. ^ formally 
diverges. This Coulomb potential is, however, screened 
by the bound valence holes and other charged impurities 
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in the system, with a screening length of the order of the 
typical Mn-Mn distance. Therefore we shall compute the 
sum in Eq. ([6]) with an exponential cut-off, AE{Rij) 
AE{R^j) e-«-/^s^, Rsc = 2 rM„, with (47rr3^„/3 = 
UMn the Mn concentration). 

Also, our calculation for the hopping matrix elements 
makes sense only for nearest neighbor sites. Correspond- 
ingly, in the Hopping part of the Hamiltonian we in- 
troduce a rigid cut-off, Rq = 2 run, to keep only hop- 
ping to the first 'shell' of atoms. Our results do not 
depend too much on these cut-off parameters; Changing 
the Coulomb energy cut-off results in an overall shift of 
the impurity and valence band energies, and does not in- 
fluence the optical spectrum. Similarly, the results are 
almost independent of Rq as long as it is in the range of, 

Rq = 2 ryin- 
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FIG. 5: Momentum matrix element [in units ft = 1] as a 
function of the separation R between the two Mn sites. For 
R > 8A the matrix elements are very well-approximated by 
the Peierls substitution (dashed line). 



2. Valence band 

In the tight-binding approach of the previous subsec- 
tion only the bound hole states appear. However, to 
describe optical transitions, we also need to account for 
extended states within the valence band. These states are 
involved in local optical transitions, where a hole local- 
ized at site TO absorbs a relatively high energy photon 
{E ~ 200meV) to make a transition relatively deep into 
the valence band. 

To account quantitatively for these transitions, we 
make the following crucial observations: (i) The tran- 
sitions involve hole states of relatively high energy, and 
correspondingly, states of a relatively short wavelength, 
(ii) Since the initial hole states are localized at the Mn 
site, the created valence band hole will also be localized 
close to it. There the major effect of neighboring Mn 
sites is to create a Coulomb potential, which is smooth 
at the scale of the wavelength of the created hole. As a 
result, we can describe the final hole state in terms of the 
following simplified 'semiclassical' Hamiltonian: 



val 2 



1 



(7) 



where — |r — Rj| is the electron coordinate measured 
from impurity i, and the shift Ai?*^[ is given by 



Ai?i val — — 



1 



^ e R,. 



(8) 



Notice that AE'i^vai ~ A£'i since most of shift AEi comes 
from the Coulomb shift in Eq. ^ (see Fig. [5]). However, 
these two energies are not exactly the same, AiJ^ vai 7^ 
AEi, since AEi also contains the kinetic shift AE, shown 
in the inset in Fig. [4l It is precisely this kinetic shift, 
which is responsible for the red-shift of the mid-infrared 
optical peak with increasing Mn concentrations, x. 

Furthermore, we observe that local optical transitions 
only involve p-states. These states vanish at the origin. 



and do not really feel the central cell correction. There- 
fore, to describe the local optical transitions, it is reason- 
able to include locally these p-states only, and describe 
the valence band in a second quantized formalism as 



band _ 
val ^ / y 



dk ( fc2 



A£;(Rj 



where a 



t 

krncT.i 



(9) 

creates a scattering state in the p-channel 
around impurity site i in angular momentum channel 
TO = x,y,x, with spin a {a = — 3/2, .., 3/2), and radial 
momentum k. The operators aj^.^^^ ■ are normalized to 
satisfy the anticommutation relation 

Wkma.i^'^k'm'fS.j} = TrSapSmrn' Sij 6{k - k') . (10) 

As mentioned above, scattering states in the p-channel do 
not feel the central cell correction, therefore the al^^^^^'s 
create to a very good approximation Coulomb scattering 
states. Notice that in Eq. ^ an independent local band 
is associated with each impurity site, i.e., propagation 
between various impurity sites within the valence band 
is ignored. This is a reasonable approximation for the 
high-energy, fast processes considered in this work. 



B. Coupling to the electromagnetic field 

At the microscopic level, the vector potential couples 
directly to the momenta of the valence holes. However, 
within the impurity band, these valence holes reside on 
specific orbitals localized at the Mn sites, and therefore, 
as a first step, we need to determine how the vector po- 
tential A couples to these impurity band states. To this 
end, we determined the matrix elements of the momen- 
tum operator p between the lowest lying states of our 
Mn2 molecule (see Appendix IB|) . We find that the effec- 
tive coupling of a homogeneous vector potential to the 
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FIG. 6: Wave vector dependence of the on-site optical matrix 
element computed for free and Coulomb scattering states [h = 
1] . We also compare this with the result obtained for a simple 
Hydrogen-like variational wave function and free scattering 
states. 



impurity band is given by 



TTimp 
-"cxt 



E 



\p{Rij)\-^A ■ nij{ic\^Cja + h.c. 



TTIqC 



(11) 



where the unit vector riy — (R^ — Rj) / Rij specifics the 
direction of the bond, and p{R) is the momentum matrix 
element extracted from the variational calculation. This 
matrix element is plotted in Fig. O As we also verified, 
for R > 8A, this tediously-obtained matrix element is 
quite well-approximated by the simple Peierls substitu- 
tion. 



p{R) w -Rt{R) TOo/7 



(12) 



The term (1111) describes intrahand transitions between 
states inside the impurity band. It is this term that gen- 
erates the Drude peak and which is responsible for the 
AC conductance. 

However, we also need to account for interband transi- 
tions, i.e. transitions from the impurity band to extended 
states in the valence band. In our approach, these tran- 
sitions are local, and are described by the Hamiltonian, 



rrtr 



E 



dk 



+ h.c.) 



(13) 

Here Am denote the m — x,y, z components of the 
vector potential A, and the on-site optical matrix ele- 
ments |p(fc)| can be computed from the variational wave 
function of a single Mn site, We determined \p{k)\ for 
Coulomb scattering states as well as for free electron 
states (details of this calculation are presented in Ap- 
pendix El ) ■ The final results are shown in Fig. [6l 

For free-electron scattering states the single-parameter 
matrix element takes on a particularly simple form, 



limn 



25/2 ^5/2fc2 

vs (a' + k^y 



(14) 



where a = 0.09lA~^. The momentum matrix elements 
vanish for fc — > since p-states have a node at the origin, 
but they also vanish for very large momenta, where the 
valence hole wave functions oscillate much faster than 
the characteristic scale of the bound state. This results 
in a maximum optical transition rate for valence states 
at about 100 meV below the valence band edge. It is this 
momentum-dependence of the optical matrix element 
that is ultimately responsible for the mid-infrared peak 
in the optical spectrum at frequencies huj « 200 meV. 
The Coulomb and free electron matrix elements behave 
qualitatively the same way and are in both cases well 
approximated by the expression obtained for a single 
hydrogen-like variational wave function, \^!). However, 
the more realistic Coulomb matrix elements have a some- 
what smaller amplitude and the maximum transition rate 
occurs at slightly higher values. In the rest of this paper 
we shall use these more realistic Coulomb matrix ele- 
ments to compute the optical conductivity. 



III. OPTICAL CONDUCTIVITY AT T = 0: 
MEAN FIELD APPROXIMATION 

In the rest of the paper we shall focus on the calcu- 
lation of the optical properties of Gai_xMnxAs at low 
temperatures. First, we shall discuss the case of T = 
temperature and treat the Mn spins within the mean field 
approximation, justified by the relatively large value of 
the Mn spins, Sun — 5/2. Later in section |V] we shall 
discuss how one can go beyond this approximation and 
compute spin wave corrections. These corrections, how- 
ever, turn out to be small, and the mean field description 
turns out to be quite accurate. This is related to the fact 
that, in the concentration range studied here, the cou- 
pling G is much smaller than the binding energy of the 
holes G ^ Eq as well as the width of the impurity band, 
and therefore it hardly influences the optical spectrum at 
these energies. 



A. Mean field Hamiltonian 

In our model Hamiltonian, Mn spins couple directly 
only to states in the impurity band. Therefore, as a first 
step, we only need to treat the coupled Mn spin-impurity 
band system, described by Eq. At the mean field 

level, the interaction part is rewritten as 

G S, cIFq = G [(S,)c1Fc, + S,(c|Fc,) 

- (S,:)(c|Fq)] + Hfinctii); , (15) 
i?fluct(i) = G (S, - (Si))(4Fc, - (4Fc,)) , (16) 

and the fluctuating part, i/fluct is neglected. With this 
approximation, Eq. Q becomes exactly solvable at any 
temperature. For given values of (Si), the hole part of 
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the Hamiltonian is diagonalized by the unitary transfor- 
mation 

i,a ^ 

where wave functions 0*(iQ;) can be found by solving 
a relatively simple eigenvalue equation. The mean field 
Hamiltonian can then be rewritten in terms of the oper- 
ators ajj and the corresponding eigenenergies, i?^, as 

^i^J = E^/^«K-Es^h,+^(s,)h,, (18) 

/i i i 

where the local fields — — G(cjFci) can be expressed 
in terms of the new basis as 

h, = -GY,F^.,{^) f{E^) , (19) 
F^^it) = ^0;(*a)F„^0,(z/3) , (20) 

a/3 

with / the Fermi function. Trivially, determines 
the expectation value (Si), that enters the eigenvalue 
equation of (j),y(ia). This field must therefore be deter- 
mined selfconsistently. At T = temperature the self- 
consistency equations become rather simple, since then 
the spins are completely polarized in the direction h^, 
(Si) — S'Mnhi/|hi|, and the Fermi function simply be- 
comes a step function. However, for small concentra- 
tions and/or finite temperature the Mn spins are not fully 
polarized, and the numerical solution of the mean field 
equations becomes time-consuming. 

B. Optical conductivity 

The mean field solution of the Hamiltonian provides us 
the equilibrium state of the coupled spin-impurity band 
system. To compute the optical spectrum, however, we 



also need to transform the coupling to the electromag- 
netic field to the 'canonical' basis, Eq. ([T7]). 

We can express the coupling of the impurity band to 
a time-dependent vector potential A(i) in terms of the 

canonical operators, as H"^^ = J^ij.^-^ ' '^jiifj^'^'^i^j 
where the current operator's matrix element is defined 
as 

[i 0*(ia)0,j(ja) - i (j>l{ja)(j)f,{ia)] . (21) 

Similarly, we rewrite the term (jl3p in the Hamiltonian in 
this canonical basis as 

^^ext= E ^™(0(jfc,m,.aL„,a, + h.c.) , (22) 

with Jk,ai,n = b(fc)| ^i(l)n{ia). 

The optical conductivity (7(w) can be computed by 
means of the Kubo formula, which expresses it in terms 
of the retarded current-current correlation function, 

a(a;) = -ilmH-}(c.). (23) 

Since, corresponding to Eqs. (|2ip and Eq. (|22p . our cur- 
rent operator consists of an interband and an intraband 
part, the current-current correlation function can also be 
expressed in terms of two contributions: an intraband 
contribution, CTintra, corresponding to transitions within 
the impurity band, and an interband contribution, CTintcr, 
describing transitions from the impurity band to the va- 
lence band. In our approach the low frequency behav- 
ior arises from impurity band scattering alone, while the 
valence band plays practically no role there. The high 
frequency conductivity, on the other hand, is typically 
dominated by interband transitions to the valence band. 
In the present, impurity band approach, these transitions 
give rise to a mid-infrared peak close to 2000cto~^. 



I 

The previously mentioned two contributions can be easily computed within the mean field approach, using the 
diagrammatic formalism presented in Section |Vl and the interband contribution can be simply expressed as 

Notice that the transition being local, it is also the local energy of a valence band hole, that appears in the denominator 
of this expression, efe(i) = + AEi. 

The intraband contribution can be written in a similar way as 

<^^ntra (c) = --^ ^-.^0^ Im | j , (25) 

with Gq the polarization of the external light. 



Clearly, both the energy and the precise structure (i.e., localized or extended character) of the impurity band 
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states enter the matrix elements in the above expressions. 
However, apart from the value of the DC conductance, 
which clearly must vanish as one approaches the local- 
ization transition, the gross high-frequency features will 
turn out to not be very sensitive to the localization tran- 
sition itself. This is not very surprising, since the lo- 
calization transition involves mostly states at the Fermi 
energy, while optical transitions are dominated by states 
deep below or high above the Fermi energy. 



IV. MEAN FIELD RESULTS 

A. Density of states 

Let us start first by discussing some details of the nu- 
merical solution of the mean field equations described in 
the previous section, and the structure of the hole states 
obtained. To start the numerical calculation, we first gen- 
erate an initial distribution for the substitutional Mn im- 
purities on an fee lattice with lattice constant a « 5.6 A. 
Nearest neighbor Mn sites are, however, not favored dur- 
ing the growth process, since Mn ions act as charged im- 
purities. To simulate this effect, we let the Mn ions relax 
through a simple classical Monte Carlo diffusion process 
while assuming a screened Coulomb interactions between 
them^ Typical values of this Monte Carlo time used are 
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FIG. 7: Ground state density of states for different Mn con- 
centrations and a fixed hole fraction, / = 0.5. The solid lines 
represent the DOS for the impurity band while dotted lines 
denote the valence band DOS. The Monte Carlo relaxation 
time is fixed to Imc = 1 per Mn ion. Shaded areas indi- 
cate locahzed states, while unshaded regions correspond to 
extended states. For these concentrations we find two mobil- 
ity edges in the impurity band: one separating states in the 
tail of the impurity band, while the other separating localized 
states in the impurity band - valence band gap. 



in the range of ^mc ~ 1 step/atom. In all our calcu- 
lations we use periodic boundary conditions to suppress 
surface effects. 

Once the configuration of the ions is determined, we 
can construct the effective Hamiltonians Eq. ([2]) and 
^ as described in Section Iff Al Having constructed the 
Hamiltonian, we then solve the mean field equations it- 
eratively and determine the equilibrium spin configura- 
tion. Here we only focus on T = temperature, where 
(S,;) = STJi, with Jli acting simply as a classical vari- 
able. In the end of the mean field self-consistency loop 
the states in the impurity and valence band are fully 
characterized, i.e. the corresponding eigen-energies 
and £k{i) as well as the wave functions are available. We 
can then compute the necessary matrix elements of the 
current operator using Eq. pdj) , and apply ([23]) and (|38l) 
immediately to compute the contributions to the optical 
conductivity. 

For a good accuracy, we consider system sizes L as 
large as possible, and we average over many configura- 
tions. Usually for a given Mn concentration and hole 
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FIG. 8: Ground state density of states and participation 
ratios for different Mn concentrations but a fixed hole fraction, 
/ = 0.5. The Monte Carlo relaxation time is fixed to Imc ~ 1- 
The shaded regions indicate the position of the mobility edge 
for each concentration. 



9 



fraction, we average over 100 different impurity config- 
urations, and we consider systems with the number of 
impurities in the range of 200. 

In Fig. [7] we present resuhs for the impurity band den- 
sity of states (DOS) for different Mn concentrations but 
for a fixed hole fraction, / = 0.5. We also indicate the 
average DOS for the valence band. Shaded regions, de- 
noting localized states were determined by by analyzing 
the finite size scaling of the participation ratio 



PR(/i) 



(26) 



For extended states PR scales with the system size, while 
for the localized states remains finite for L — ^ oo. Fig. [8] 
shows the details of such a finite size scaling analysis. 

As we expect, increasing the Mn concentration, the 
impurity band gets broader and starts to overlap with 
the valence band for Mn concentration as low as 0.5% 
(see also Figl?]). However, at the same time increasing 
Coulomb disorder shifts the tail of the impurity band 
deeper inside the gap (the top of the valence band is 
fixed at Ey^ = 0). States inside this tail are strongly 
localized and have an overwhelming localized impurity 
state character. 

We clearly observe a metal insulator transition (MIT) 
at a critical concentration x « 0.2%, which corresponds 
to a hole concentration p ~ 4.5 x lO^^cm"'^ for / « 1. In 
our approach this metal-insulator transition is simply an 
Anderson localization transition, since electron-electron 
interaction effects are neglected: For x < 0.2% all the 
states in the impurity band are localized by disorder and 
the system is an insulator, while for x > 0.2% the metal- 
licity of the system still depends on the position of the 
Fermi level inside the impurity band. As an example (see 
FigUl bottom layer) for x = 0.3% and hole fractions less 
than / = 0.8 the system is still an insulator, while for a 
hole fraction / ^ 1.0, the Fermi level is already inside the 
delocalized region and we find a metallic state. Increas- 
ing the Mn concentration to larger values, the extended 
region grows and the system has a metallic behavior for 
typical values of the hole fraction, / ~ 0.5 1, (see also 
Fig. [T]). For hole fractions / « 0.6, e.g., we find that the 
ground state is metallic for x > 0.5%. These numbers are 
in very good agreement with experimental results, where 
the MIT has been reported to occur in a concentration 
range x « 0.2 - 0.3%M^ 

We must remark that here we completely neglected 
electron-electron interaction. Electron-electron interac- 
tion can obviously modify many of the physical properties 
close to the MIT and lead to the appearance of a Coulomb 
gap as well as Altshuler-Aronov-type anomalies. ■^^ On the 
other hand, even annealed Gai_xMnxAs samples are typ- 
ically compensated, and the electron-electron interaction 
has been found to be less relevant for such systems. Also, 
we can argue that our results are in a way self-consistent, 
and probably not very sensitive to Coulomb interaction; 
For typical concentrations x < 3% and hole fractions in 
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FIG. 9: Shifted Local density of states for different Mn con- 
centrations. The Monte Carlo relaxation time is fixed to 

tMC = 1- 



the range / ~ 0.3-^1, we find that aproximately 60-;-80 % 
of the sites have less then one electrons on them, and 
less then ~ 5 % of sites are more then doubly occu- 
pied. Furthermore, even for the most localized orbitals, 
the participation ratio is rather large ^ 6 — 10. There- 
fore, a typical site is occupied by a single electron and 
the Coulomb correlations are not too relevant. Treat- 
ing the Coulomb interaction appropriately is, however, 
a real theoretical challenge even for much simpler model 
systems, and is beyond the scope of the present work. 

Let us close this subsection by presenting a rather cu- 
rious quantity that we termed 'shifted local density of 
states' (shiftcd-LDOS), where at every site we computed 
the local density of states by subtracting the binding en- 
ergy of a hole, Eq, as well as the value of the average 
Coulomb shift. While the unshifted LDOS is rather fea- 
tureless, this quantity displays a sharp peak around zero 
energy. A more detailed analysis reveals that this sharp 
peak is associated with localized states having a small 
participation ratio, i.e. states in the tail of the impurity 
band. To show this, we also plotted the contribution of 
states with PR's smaller than a given value to the shifted- 
DOS. Clearly, the peak observed at zero energy is entirely 
due to localized states deep inside the gap. The physical 
interpretation of this peak is thus plausible; States in the 
tail of the impurity band are generated by random and 
large fluctuations of the Coulomb disorder. These states 
are simply localized impurity states deep inside the gap 
that are being mostly occupied, and they give rise to lo- 
cal optical transitions which have a strong impurity state 
transition character. 



B. Optical conductivity 

Having diagonalized the mean field Hamiltonian, it 
is straightforward to compute the optical conductivity. 
As we shall see, our theory, which incorporates inter- 



10 




Frequency (cm ) 

FIG. 10: The real part of the optical conductivity for different 
Mn concentrations and hole fractions. In a) and c) the Mn 
concentrations is 0.03%, while in b) and d) it is a; = 1%. 
In c) and d) the intraband and interband contributions are 
presented for a hole fraction / = 0.5. The sample in panel 
a) is insulating while the one in panel b) is metallic. The 
peak observed at ~ 2000 -ir 2500cm~^ for the metallic sample 
is due to the valence to impurity band transitions while the 
broad feature at smaller energies is a reminiscent Drude peak 
generated by the impurity band contributions. 



impurity band transitions as well as intraband transi- 
tions, accounts qualitatively as well as quantitatively for 
all major features of the experimentally observed opti- 
cal spectrum. We shall also present a scaling analysis 
of the optical conductivity close to the metal-insulator 
transition and compute the dynamical critical exponent. 

The Drude contribution, i.e. the lu feature is 
controlled by excitations close to the Fermi level. Since 
in this small concentration limit the Fermi level resides 
in the impurity band, the "Drude peak" and also the 
DC conductivity are expected to be dominated by the 
intra impurity band contribution. The contributions of 
extended valence band states is negligible for a; <C i?o, 
since these states are too far from the Fermi level to be 
of any relevance for the Drude contribution. 

The interband contribution, on the other hand, is gen- 
erated by transitions between the impurity and valence 
bands. As can be seen in Fig. [H the two bands start to 
overlap for Mn concentrations x > 0.5%. However, al- 
though it also depends on the hole fraction /, the energy 
gap between the Fermi level and the top of the valence 
band is always in the range A = O.le^. We therefore 
expect to see a broad mid-infrared feature in the energy 
range A < ui < 2A. 

These expectations are indeed met by our numerical 
results shown in Fig. [TOl In panel d) we also display the 
intraband and interband contributions separately for a 
typical Gai-xMuxAs sample in the metallic regime, with 



X = 1% and a hole fraction / = 0.5. The intraband con- 
tribution has a Drude-like behavior, while the interband 
signal presents a peak at the energy 2500cm~^. How- 
ever, the Drude contribution is rather broad, and in fact, 
already for this relatively small concentration, it almost 
completely merges with the interband contribution. The 
overall behavior is in good agreement with the experi- 
mental data presented in Ref.[9 10 11 12]. Remarkably, 
not only the peak position but also the overall magnitude 
work out to have values in the experimental range. 

As already shown in Fig. [31 increasing the number 
of carriers leads to a red-shift of the mid-infrared peak 
in agreement with the experimental observations; more 
metallic samples tend to have a strongly red-shifted mid- 
infrared peak at a frequency, ~ 1500 cm~^ rather than 
at ~ 2500 cm~^. We find that in the concentration range 
considered here, to a good approximation, the red shift 
only depends on the effective optical spectral weight, 
iVefi, defined through Eq. (|28p. Remarkably, there is 
an almost perfect match between our theoretical results 
and the experimental values for intermediate concentra- 
tions. The decreasing trend observed in the position of 
the mid-infrared peak as function of hole fraction can be 
understood simply from Fig. [3l for a given Mn concen- 
tration, increasing the carrier numbers shifts the Fermi 
level closer to the top of the valence band. As a conse- 
quence, the optical gap decreases and the mid-infrared 
peak moves towards lower energies. 

A first look at the results presented in panels c) and d) 
of Figure [10] can drive us to the conclusion that there 
is not much difference between the overall optical re- 
sponse of the metallic and insulating samples. However, 
the effective mass analysis and spectral weight analysis 
presented in the following subsection clearly shows the 
difference between these phases. 

However, before turning to the effective mass analysis, 
let us discuss the low frequency properties of the opti- 
cal response close to the MIT transition. The scaling 
theory of Abrahams et alM- has been extended to the 
dynamical conductivity by Shapiro et al.^ who showed 
that at the mobility edge the ac conductivity of a d = 3- 



dimensional system obeys a power law, <t{uj) 



They also found that on the metallic side the conductiv- 
ity goes as — (7(0) ~ ujid.-'^)/'^ while in the insulating 
phase (j{uj) ~ uj"^. 

To our knowledge, the first numerical evaluation of the 
optical conductivity close to the localization transition 
based on the Kubo formula was done in Ref . [i^l , where 
indeed it was shown that the dynamical exponent for 
the optical conductivity of a three-dimensional Anderson 
model at the critical point is 1/3. More recently, the anal- 
ysis was extended to unitary and symplectic systems^ 
where the same exponent was found. In the lower panel 
of Fig. [11] we present the low frequency limit of the intra- 
band optical conductivity close to the MIT transition for 
a fixed hole fraction / = 0.8, where the critical carrier 
concentration is found to be Pc — 4 x lO^^cm"'^ [x — 0.2) 
by a participation ratio analysis. For p > pc the low 
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FIG. 11: Upper panels: a) the dc-limit of the intraband 
optical conductivity for concentrations between 0.02% and 
1%. b) The scaling exponent a for the conductivity in the 
low energy limit as function of the carrier concentration. The 
MIT transition is indicated by a drop of the critical coefficient 
a down to 0.35. Bottom panel: Frequency dependence of the 
intraband optical conductivity in the small frequency limit 
presented for different Mn concentrations. In all three figures 
the hole fraction was fixed to f = 0.8. 



energy optical data can be nicely described by a power 
law, (j{lu) = a{0) + Au!°'. Since we use periodic boundary 
conditions, we observe a continuous metal- insulator tran- 
sition similar to the case of the Anderson model42. For 
an infinite system, the dc conductivity cr(0) approaches 
zero at the transition point and vanishes in the insulating 
phase. Of course, for a system of a finite size like ours, 
a{0) is not strictly zero even in the insulating phase, and 
only a crossover is observed between these two regimes 
(see the upper left panel of Fig. [TT|) . 

The most relevant quantity is the exponent a. In the 
metallic state a best fit gives a ^ 0.5 in good agreement 
with the analytical predictions. Although the error bars 
are rather large, a seems to slightly decrease down to 
0.35 at the MIT where it starts to increase again in the 
insulating regime. The two approaches, the participation 
ratio and the scaling analysis of the dynamical conduc- 
tivity based on which we have identified the MIT point 
are in good agreement. In both cases the critical point 
is associated with approximately the same hole carrier 
concentration. 



C. Effective mass analysis 

Experimentally, optical conductivity data are often 
used to extract important information on the charge car- 
riers such as their concentration or effective mass, or the 



value of kpl. From the analysis of the experimental data, 
two types of effective masses can be extracted: (i) the ef- 
fective mass m* that enters the Drude formula for the 
resistivity, i.e. the dc limit of the optical conductivity. 
(m) the optical mass Wopt that is related to the spectral 
sum rule. These two masses are often quite different. 
Also, they are typically extracted based upon some as- 
sumptions (isotropical mass and quadratic Hamiltonian, 
etc.), which are violated in the real system. Nevertheless, 
they both contain important information, and it is there- 
fore interesting to determine them from our theoretical 
results and compare to the experimental data. 

The effective mass m* can only be roughly estimated 
from the data. Following the procedure of Ref [ll|, we 
typically find it to be in the rage of m* 10 — 30 mg. 
The mass m* is also related to the the metallicity param- 
eter, kpl, with kp the Fermi momentum and I the mean 
free path. In fact, this dimensionless parameter can be 
extracted much more reliably from the experimental data 
than m* , and it has a much more obvious interpretation: 
Small values of kpl correspond to samples in the insulat- 
ing phase or close to the metal insulator transition, while 
good metals are characterized by large values of kpl. We 
shall therefore focus on this parameter. 

The value oi kpl can be estimated from the hole con- 
centration and the dc conductivity using the Drude for- 
mula, rewritten as 



kpl^ i ^ 



a(0)^A. 



(27) 



Here s is the spin degeneracy of the band. In an s = 4 
four band model we found that for carrier concentrations 
between p — 2.0x 10^" 9.0 x lO^'^cm^^ our optical con- 
ductivity would correspond to values 0.3 < kpl < 3.0. In 
other words, for these high resistances the mean free path 
is already less than the Fermi wavelength Xp = 2iT/kp. 
Thus the Drude analysis clearly leads to a inconsistency, 
and shows that the optical conductivity does not origi- 
nate from a weakly perturbed valence band. While in the 
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FIG. 12: The optical mass as function of number of carries. 
Only the impurity band contribution was used to compute 
niopt- Deep in the metallic regime, for carrier concentrations 
larger than 2x lO^^cm"^ the effective mass remains practically 
constant, while in the critical regime as approaching the MIT 
transition (p ~ 10^^ -r- 10^" cm'^) an power- like increase in the 
effective mass is observed. 
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present work, where we use an impurity band approach, 
this is obvious, the same analysis can also be performed 
for the experimental data to show, that a weakly per- 
turbed valence band picture is inappropriate to describe 
most Gai-xMuxAs samples.^ The dc conductivity values 
obtained within our impurity band picture are, on the 
other hand, perfectly consistent with the experimentally 
extracted values of kpl. 

The optical mass, Wopt, can be defined through the 
effective optical spectral weight, N^g, defined by the in- 
tegral 

p 2 f^" 

= iVoff = 5" / CTintra (t^) du , (28) 

with ujc a somewhat arbitrary energy cut-off, which has 
been set to 800 meV in Ref. [11]. The optical mass rriopt 
defined through this formula can be quite different from 
the microscopic mass of the carriers, mo, and it provides 
a useful measure of the "heavyness" of the charge carri- 
ers. To compare with the experiments, we extracted rriopt 
from our data, and in Fig. [12] we plotted it as a function 
of carrier concentration. Deep in the metallic regime, 
we have obtained an optical mass that remains approxi- 
mately constant, rriopt /rrie ~ 1, for carrier concentrations 
larger than p « 2 x 10^", i.e. for Mn concentrations larger 
than about x ~ 1%. Our kpl values as well as our optical 
mass results are in good agreement with the experimen- 
tal data,— where the optical mass was found to be in the 
range 0.7 < rriopt/me < 1.4. Notice that these values 
are about twice as large as the bare valence band mass, 
mo = 0.56 rrie. As we approach the MIT transition, 
mopt increases rapidly, and can reach values as large as 
rriopt ~ 10 rrie. This large optical mass renormalization 
is characteristic of the vicinity of the critical concentra- 
tion, (see Fig fT2)) . In this regime we find a metallicity 
parameter kpl that is always less than 0.5. 
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TABLE L Vertices functions within the interacting field the- 
ory. Dashed lines represent impurity hole propagators, while 
wavy lines denote bosonic propagators. 
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(29) 



where we have neglected terms that are subleading in 
1/5. 

In this language, the mean field Hamiltonian, Eq. JT^ 
simply becomes 



Hl^l = Y,E,,ala^ + J2\i^Mk, (30) 
V. SPIN FLUCTUATIONS u 



So far we have treated Mn spins at the mean field level. 
In the present section we shall discuss how one can go be- 
yond mean field by systematically including spin fluctua- 
tions. To do that, we shall make a 1/ S expansion around 
the classical limit, and represent spin fluctuations using 
Holstein-Primakoff bosonsi^ The ground state is, how- 
ever, generally non-coUinear j^^i"^^ leading to some com- 
putational complications. 

Here we shall consider again only the T ~ temper- 
ature limit. There the mean field equations imply that 
(Si) = ^Mnef, with the unit vector ef = hi/|hi| being 
parallel to the external field created by the valence holes. 
We shall use this direction as a quantization axis of the 
Mn spin at site r, and introduce two additional unit vec- 
tors, e^'^ perpendicular to e^. Using Holstein-Primakoff 
bosons, we can then represent the Mn spin operators at 



while corrections to the mean field the arise from the 
fiuctuation part of the interaction, Eq. pB]) , neglected so 
far: 

i?fluct(r)«-G6t5, ictefFc, : (31) 
+ G ^{bl-.cle+Far.+h-.clerFc. :) , 

where we defined = ef ± ief , and : • • • : denotes 
normal ordering with respect to the mean field ground 
state. Thus the neglected fluctuation terms just describe 
the interaction of spin waves with spin excitations of the 
impurity band. 

We can also rewrite the interaction part of the Hamil- 
tonian using the proper eigenstates of the mean field 
Hamiltonian and the corresponding creation and anni- 
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hilation operators as 

+ (^^(*) ^l^'^^l + ala^b, 

with the couphngs 7™ defined from Eq. (|^D)) as 



(32) 



G 



G 



(33) 



Having introduced a bosonic representation for the 
spins, we can now use standard diagrammatic methods 
to compute spin-fluctation corrections to the optical con- 
ductivity perturbatively in the small coupling, G. In the 
noninteracting theory, defined by Eqs. (O and pO)l . we 
have three fields, and, correspondingly, we have three 
Green's functions. The non-interacting Green's functions 



lUJn — i^u 



Gta|(fc,J,^n) = - 



1 



iujn - £k{i) 



(34) 
(35) 



describe single particle excitations in the impurity and 
valence bands. We shall denote them by dashed and 
continuous lines, respectively. Spin waves are described 
by the noninteracting bosonic propagator, 



1 



(36) 



This propagator is site-diagonal, and accounts for the 
spin precession created by the local exchange field of the 
impurity band holes. The interaction part, Eq. (|32p . 
introduces then three types of vertices between these 
Green's functions, which we depicted in Table [J 

To compute the current-current response function, we 
shall use a random phase approximation-type (RPA- 
type) approximation in the bosonic line, which mediates 
the interaction between the valence band holes: In the 
diagrammatic language this means that one neglects self- 
energy corrections, and sums up only the diagrams shown 
in Fig. [T3] (top layer). Notice that at T = the vertex 7^ 
does not give a contribution to the series in leading order 
in G. The RPA series can be converted into an integral 
equation (Bethe-Salpeter equation), as shown in Fig. 1131 
(bottom layer), where one formally defines a renormal- 
ized current vertex function, F^^. Expressed analytically, 
the vertex equation reads 



E 



^7;.c(*a.,H.)n^,(zao 



^ fj.v — j^iv (37) 




FIG. 13: Top: RPA series giving the leading correction to the 
optical conductivity. The dot represents the current vertex. 
Middle: (a) interband contribution to the optical conductiv- 
ity and characterize excitations, (b) intraband contribution. 
Bottom: Bethe-Salpeter equation for the renormalized cur- 
rent vertex P^i/, indicated as a filled .triangle. 



Here P^,y(ir2„) stands for r^y(ii7„) = r^y(ir2„)eo with 
Go the polarization of the electric field, and similarly, 
— i^iv^o- The intraband polarization bubble is then 
expressed in terms of the bare j^j/ and full vertex P^^ as 

T^TriiVLr,) - JM-r^^(zaO no_,(zf^„),, (38) 

with the bare polarization buble defined as: 

KA'^n) = E ^^n) G^^{v, zc.„+zJ7„) (39) 

Immediately, this gives for the intraband optical conduc- 
tivity the expression 

a^ntra {^) ^ 0-^60) (P,,. (f^)eo ) [f {E ^) - f{E,)] 



n + E^,-E^ + i5 VI- E^ + E^ + i5 



(40) 



The mean field limit is recovered if one sets 7 ^ in 
the Bethe-Salpeter equations, Eq. ([27]). Then Eq. ([iD]) 
reduces to Eq. (P5)l . 

We have solved the Bethe-Salpeter equations Eq. ([57]) 
numerically, and computed the resulting corrections to 
the optical conductivity. However, we found that, at 
typical optical frequencies, they do not give an impor- 
tant correction to the optical conductivity as compared 
to the previously presented mean field results. This is 
not very surprising: Magnetic excitations have energies 
in the range of Tp, i.e. in the 5-6 meV range, which is 
very small compared to the optical frequencies studied 
here. In the frequency range oj < 5meV they are, how- 
ever, expected to result in interesting effects, and they 
may lead to some additional features in the optical spec- 
trum. 
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The present calculation has been carried out at T = 
temperature. In this case the diagrammatic derivation 
of Eqs. ([37| and ((40|) is very transparent and straight- 
forward. However, we also derived these equations using 
the much more tedious equation of motion method of 
Ref. [45], which holds also at finite temperatures. It 
turns out that Eqs. (|37p and (I40p are very general, and 
in the same form, they carry over to finite temperatures 
too. The only modification is that the original mean field 
Hamiltonian must be diagonalized using finite tempera- 
ture expectation values. Our formalism thus provides a 
way to study the effects of finite temperature on the op- 
tical conductivity and the interplay of ferromagnetism. 
The study of these finite temperature small frequency 
effects is, however, beyond the scope of the present work. 



VI. CONCLUSIONS 

In this paper, we presented a calculation of the opti- 
cal properties of Gai_xMnxAs in the very dilute limit, 
where it can be described in terms of an impurity band 
picture. Our approach consisted of constructing an effec- 
tive Hamiltonian with parameters determined from mi- 
croscopic variational calculations. The effective Hamilto- 
nian obtained this way not only accounts for the impurity 
band, but it also describes transitions from the impurity 
band to the valence band. 

Our mixed approach captures correctly the most essen- 
tial features of the experiments: It predicts a rather wide 
Drude peak, originating from inter-impurity band tran- 
sitions, which is partly merged with a mid-infrared peak 
at ~ 200 meV. The overall conductivity values as well as 
the positions of these features are well-reproduced by our 
theory. Remarkably, our calculations give a quantitative 
description of the concentration-induced red-shift of the 
mid-infrared peak, and gives optical mass and residual 
resistivity values that agree well with the experimentally 
observed values. 

Furthermore, we are also able to capture the metal- 
insulator transition, which, depending on the level of 
compensation, occurs in the range of a; = 0.1 — 0.3 % 
for uncompensated samples while it is in the range of 
X — 1 ^ 2 % for moderate compensations. This is in 
good agreement with the experiments. Our results thus 
promote the theoretical picture that the metal-insulator 
transition takes place within the impurity band, well be- 
fore it merges with the valence bandi^i^ This picture 
is indeed supported by the observation of a finite ac- 
tivation energy as one approaches the metal-insulator 
transitioni^i Further indirect evidence for this scenario is 
given by Fig. 8 of Ref. [l^: at very small concentrations, 
where the impurity band is expected to be insulating, 
and the conduction is due to activated behavior to the va- 
lence band, the mobility of holes agrees with that of many 
other alloys, that are known to have valence band conduc- 
tion. However, right after the metal-insulator transition 
(x > xc ~ 0.2%), the mobility clearly drops to a much 



-o- 




Valence band 



FIG. 14: Deep, strongly localized occupied levels in the tail of 
the valence band move together with the valence band edge, 
and can give rise to impurity transitions at about 200 meV, 
indicated by the wavy line. Empty circles denote holes. 



smaller value. The strong deviation from other 'valence 
band' compounds could naturally be explained by the 
fact that conduction is due to holes within the impurity 
band, and that the metal-insulator transition also occurs 
there, in agreement with our calculations. 

We have also investigated the scaling of the conduc- 
tivity and the optical mass close to the metal-insulator 
transition. We find that the intraband optical conduc- 
tivity scales as o'(w) = <t{0) + Auj^^^ on the metallic side, 
close to the transition. Close to the transition point 
the optical mass is strongly renormalized and takes val- 
ues mopt ~ 10 me, while far from the metal- insulator 
transition we find a mass in the range Wopt ~ rUe- 

Although our calculations were performed at the mean 
field level, we extended it to incorporate spin wave fluc- 
tuations too. While the formalism discussed here is only 
applicable at T = temperature, we can also derive 
the same integral equations using an equation of motion 
method (not discussed here) that carries over to finite 
temperatures. The final equations presented here thus 
apply for finite temperatures too. Nevertheless, we find 
that for typical optical frequencies spin waves do not give 
a substantial contribution. They may, however, influence 
the low frequency optical properties (hcu < 20 meV) . 

Our approach is designed to work in the limit of small 
concentrations, and it should break down for large con- 
centrations. Where exactly this breakdown occurs, is 
not quite clear. In Ref. [iJI it was argued that the 
break-down should appear at some active Mn concen- 
tration in the range of x ~ 1%.— However, surpris- 
ingly, all spectroscopic data seem to favor a picture in 
terms of impurity band physics even at intermediate 
concentrations ^iS dOi n 1 13 Also, our results are in very 
good agreement with the experimental optical data, if 
we just blindly extrapolate them to a: ~ 3 — 4%. The- 
oretically, this is rather mysterious: One possible expla- 
nation would be that Gai_xMnxAs is rather inhomoge- 
neous, and some spectroscopic data are dominated by 
regions of small Mn concentration. The other possibility 
is that although the impurity band density of states may 
merge with the valence band at larger concentrations, 
most of the occupied states are in the tail of the valence 
band, which has a strong impurity band character. This 
scenario is sketched in Fig. [TH Since optical transitions 
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are local for these states, and since the valence band and 
the deep acceptor states are Coulomb shifted by approx- 
imately the same amount, these levels deep in the gap 
give rise to an impurity state transition at a frequency 
close to ^ 200 mcV. This is indeed supported by our cal- 
culations, where wc find that the shifted local density of 
states has a large peak coming from localized states in the 
tail of the impurity band. A third possibility would be 
that the properties of Gai_xMnxAs on its surface, tested 
by essentially all spectroscopical probes, differ from those 
of the bulk, and thus the effective concentration of Mn 
ions is somehow reduced there. 
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APPENDIX A: GROUND STATE WAVE 
FUNCTION OF A SINGLE Mn ION 

In this appendix we present some details on the varia- 
tional calculation of the bound acceptor state of the Mn 
impurities and the computation of the corresponding op- 
tical matrix elements. 



The full variational solution computed using 50 param- 
eters ai is plotted in Fig. 1151 where it is also compared 
to a single parameter variational solution. The single- 
parameter solution with a — 0.091 gives a remarkably 
accurate approximation for the wave function, excepting 
the regime r < tq, where small deviations appear due to 
the central cell correction. 
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FIG. 15: Normalized ground state wave function: Compari- 
son between the single parameter and the full calculations. 



2. Calculation of on-site optical matrix elements 



1. Acceptor wave function 

We describe the Mn acceptor by the following Hamil- 
tonian (h — nio — e — 1), 



2 er 



(Al) 



Here the factor 7 = 1.782 describes the mass renormal- 
ization term, and e = 12.65 is the dielectric function for 
GaAs. The explicit form of the central cell correction, 
Vcc{r) was given in Eq. (|4]). 

To compute the ground state wave function of (lAip . 
we used the following variational Ansatz 



Once the variational wave function (|A3p is known (at 
hand), we can compute the optical matrix elements to 
valence band states. We have two obvious choices: The 
simplest way to estimate these matrix elements is to ne- 
glect the 1/r term as well as the central cell correction 
(which does not interact too much with the p-channel), 
and to use free electron states 



(A5) 



with ji{kr) the spherical Bessel functions and Yim{0,'p) 
spherical functions. Alternatively, we can neglect only 
the central cell correction and use Coulomb scattering 
states, 



*(r) = ^A,V(a»,r) , 



1=1 

(rla) = ^a^/^e""'- 



(A2) 
(A3) 



In this expression we fixed the coefficients a's and used 
only the Ai^s as variational parameters. The variational 
equation for the latter is simply a linear equation of the 
form: 

n n 

Hc,a, Ao,^ ^EY^ Sa^a, ^a, • (A4) 

For the Hydrogenic wave functions used, the matrix ele- 
ments Haf3 = {a\H\f3) as well as the overlap parameters 
Sap = (a 1/3) can be computed analytically. 



^/2 



Fi{kr)Yirn{e,ip), 



(A6) 



where Fi is related to the confluent hypergeometric func- 
tion, ii^i(a|6|a;), as 



Fi{kr) ^ ci{r])e"'''{kry+' iFi{l + 1 + iT]\2l + 2\ - 2ikr) , 

(A7) 

with ?7 = the constant q(?7) = 2'e-"'/^ '^[f+^+^'j'^l , 
and r(z) the gamma function. Both scattering states are 
normalized to satisfy {l'mk\l'm'k') = ■nSu'Smm'^ik — k'). 

The matrix elements of Pz can be computed analyti- 
cally in both cases. For free valence electrons we obtain 



(l,0;fc|p,|vI/)frec = ^25/2VA 



V3 



i=l 



(A8) 
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while for Coulomb scattering states we have 

(1, 0; k\p,\^)coui = -^25/2e-''/2 |r(/ + 1 + 2r^)\ 

The matrix elements obtained this way were presented 
in Fig. [6] The single parameter ground state wave func- 
tion gives a very good approximation in both cases. The 
overall energy dependence of the transition matrix ele- 
ments is qualitatively similar in Coulomb and the free 
scattering state approximations. However, the height of 
the peak is about a factor of two larger in the free electron 
approximation . 

APPENDIX B: TWO SITE PROBLEM 




R [Angstrom] 



In this Appendix we determine the energies of the 
molecular orbitals of an Mn2 system and use them to 
calculate the parameters of an effective second quantized 
Hamiltonian. We consider two Mn ions, located at po- 
sitions r = Ri_2 — ±R/2, and described by the Hamil- 
tonian of Eq. ([3]). Similar to the Mn ion case, we first 
construct the lowest lying 'molecular states' using varia- 
tional wave functions of the form, 

N 

E E E^^A.,^^.^'^("^'^)' (Bi) 

fj.—s,p i— 1,2 l—l 

constructed from Hydrogen-like s and p-type wave func- 
tions centered at the Mn sites, Ri,2 

V^f ) (a, r) ^ (r |z, s, «) = -^a'/^e-'^^' , (B2) 

V-i"^ (a, r) ^ (r I*, p, a) = -l=a5/2^^g-an/2 (gg) 
4v 27r 

Here x is the direction connecting the two sites, and — 
r — Ri denote the position of the valence hole relative 
to Ri,2- We considered only states that do not have a 
mirror plane that contains x, and therefore included only 
Px orbitals. Similar to Appendix]^ we fix the parameters 
aj in Eq. (jBip . and consider only the Ai^ as variational 
parameters. 

As in case of the single Mn ion problem, the coefficients 
Ai^ ^ j satisfy a set of linear equations, 

N 

E E E[^ir'("'='"0-i?5fr^K,«o]^.,.., 

iy=s,p j=l,2 1=1 

(B4) 

with the overlap matrix elements and the Hamiltonian 
matrix elements defined as 

S'^'^j"\a,f3) = (i,^,a|j>,/?) , 



FIG. 16: The lowest few energy levels of the Mn2 molecule as 
a function of distance. The left figure presents the case where 
only s- waves were considered (18 parameters were used). The 
right panel shows the results when both s and p waves are 
considered (14 parameters were used). 

The above integrals can be performed analytically using 
elliptic coordinates. Some details of the derivation as 
well as the analytical expressions are given as an on-line 
supplement. 

We remark here that Eq. (jB4p not only provides an 
accurate estimate for the ground state energy, but it also 
accounts for the first few excited states. In Fig. [16] we 
show the first few lowest energy states obtained in this 
way as a function of Mn separation, R. As demonstrated 
in the figure, the inclusion of p-orbitals does not result 
in a major improvement for the first five states, though 
it introduces a level crossing between the first and the 
second excited states at a distance R ^ GA where the 
tight-binding approximation fails. 

The hopping t and the energy E appearing in Eq. ^ 
and shown in Fig [3] can be extracted from this spectrum 
as: 

t{R) ~ -^{Eodd{R) ~ Ecvcn{R)), 
E{R) = -{EcYcn{R) + Eodd{R)), 

where E^vcn and -Bodd represent the lowest lying (even) 
state of the Mn2 molecule and that of the first (odd) 
excited state. 

The electric field induces transitions between the 
previously-mentioned even and odd states, and thereby 
generates a coupling of the form, Eq. ifTTj) . To deter- 
mine the parameters appearing in this equation, one 
simply needs to compute the momentum matrix ele- 
ments ('(/'odd(?')|p|V'cvon(?')) that enter the coupling to 
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the external electromagnetic field. Since we aligned 
the two Mn ions along the x direction the only non- 
vanishing component of the momentum matrix elements 
is (V'odd('')|Pa; I V'cvcn (''))■ K is possiblc to cxprcss this ma- 



trix element analytically for our variational wave func- 
tions, though the final expressions are rather cumber- 
some. The final results have been plotted in Fig. [H 
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A careful calculation using the full six band Hamiltonian 
shows that the hopping is indeed fairly isotropic as long as 
the hopping distance is less than about R ~ 15A .— 
For compensated Gai_xMnxAs many of the substitutional 



Mn ions are believed to be inactive and not to participate 
in the ferromagnetic state due to the presence of interstitial 
Mn ions. 



